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Abstract 

We have considered the penguin interaction contribution to K — > 
2n decays. In particular, we have investigated the effect of the 
momentum dependence of the penguin coefficient. Our analysis 
is performed within the Chiral Quark Model where quarks are 
coupled to the pseudoscalar mesons, which means that hadronic 
matrix elements can be calculated in terms of quark loop dia- 
grams. 

We have inserted the momentum dependent penguin coeffi- 
cient into the relevant quark loop diagrams for K — > ir. We dis- 
cuss two possible prescriptions for performing the calculations, 
and conclude that the momentum dependence of the penguin 
coefficient increases the amplitude by 10-20 %. In any case, the 
(CP-conserving) penguin contribution is very sensitive to the val- 
ues chosen for the involved parameters. 



1. Introduction 



The penguin (gluonic monopole) diagram was proposed as an explanation 
of the AI = 1/2 rule |IJ. This diagram induced pure AI = 1/2 four quark 
operators into the effective Lagrangian relevant for K — > 2-k decays and 
other AS = 1 transitions. In general, non-leptonic decays are described by 
an effective Lagrangian [§] 

C(AS = 1) = ~G F K CiQ, , (1) 
V2 i 

where Gf is Fermis coupling constant and A u = V SU V^ U is the relevant 
Kobayashi-Maskawa[§] factor for CP-conserving AS = 1 transitions. In such 
effective Lagrangians, the heavy mass scales are integrated out and their ef- 
fects are contained in the Wilson coefficients Cj. The quark operators Qi 
involve the three light quarks q — u,d, s. The coefficients will in general in- 
clude short distance QCD effects calculated by means of perturbation theory 
and the renormalization group equations (RGE). Both the coefficients Cj and 
the matrix elements of the quark operators depend on the renormalization 
scale (jl (taken to be of order 1 GeV), in such a way that the physical processes 
are independent of fi. Concerning the discussion on the importance of the 
penguin interaction for the AI =1/2 rule within the standard approach, we 
refer to the literature|| f|, ||, |, |7|, |, llPf . In this paper, we will revisit this 



issue from some more unconventional point of view, and ask if some effects 
have been left out. 

To calculate the matrix elements of quark operators between physical 
hadronic states is in general a difficult task, and one normally uses various 
models or assumptions. In this paper we will use the Chiral Quark Model 



(xQM), advocated by many authors [0, |12[ [13|]. The model consists of 
the ordinary QCD Lagrangian and in addition a term C x . This new term 
includes the Goldstone meson octet in a chiral invariant way, and provides 
meson-quark couplings that makes it possible to calculate matrix elements of 
quark operators in terms of quark loop diagrams. The model is thought to 
be applicable for quark momenta below some scale of the order A x , the chiral 
symmetry breaking scale. Although this is a model, we find it very interesting 
because it reproduces the chiral Lagrangian terms for strong interactions 
to good accuracy to order p 4 (T^|. As a field-theoretical model, it can be 
formulated in terms of path integrals. Integrating out the quarks, one obtains 
chiral Lagrangians. Alternatively, integrating out the Goldstone mesons, one 



obtains Nambu type models [D 

In order to predict physical amplitudes like K —>■ 2tt, one has to match 
the high energy description, contained in the coefficients Cj, with the low 
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energy description. This matching should be done at some scale /i where 
both descriptions are valid. How to do this has been a delicate question, 
and especially for the penguin contributions because the coefficients are very 
sensitive to the quantity m c / /i, where m c is the charm-quark mass. 

In this paper we will focus on the following aspect of the effective La- 
grangian ([]]): The coefficients Ci in ([!]) originates from some Feynman dia- 
grams and will in principle depend on the external quark momenta, that is, 
the quark momenta corresponding to the quark fields in Qi in (|IJ). Working 
within the xQM, we might keep the full momentum dependence of the pen- 
guin coefficient, calculated to lowest order in the strong coupling constant. 
Being a field theoretic model, the xQM can of course account for such mo- 
mentum dependent effects. Unfortunately, our calculation will have some 
limitations. Especially, we lose the RGE analysis of the Wilson coefficient 
because the RGE analysis does not keep track of the individual quark mo- 
menta, only their overall scale. Thus, our procedure will be less systematic 
than the standard one p], |[ §, |], Still, because the standard method 
might have some difficulties with the matching around 1 GeV, we find it 
fruitful to perform such a calculation to see if there are effects which are lost 
in the standard approach. 



2. The penguin interaction at quark level 

The penguin diagram is shown in Fig.l, and induces an effective interac- 
tion (in the CP-conserving case) 

Cp = —V2Gf^ u CpQp (2) 

where Cp is the Wilson coefficient of the operator 

Q P = (d lv Lt a s) qYt a q , (3) 

where a sum over q — u, d, s is understood. The t a, s are the colour matrices 
(for a = 1, ..,8), and g s is the strong coupling. For the CP-conserving case 
considered here, the dominating part is due to u- and c-quarks running in 
the penguin loop, and the following expression is obtained in the leading 
logarithmic approximation before any RGE analysis is performed: 
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C P = -£ln^ (4) 



a s , m 
— In — 

37T /J, 

where m c is the c-quark mass. Going beyond the leading logarithmic approx- 



imation, one finds the following expression jig , |10| : 



where a s = g 2 /47r. Here the constant next to leading term is numerically 
as important as the leading logarithm. This illustrates the difficulties in 
treating the penguin diagram perturbatively. 

The expression (15]) still does not take into account the variation of Cp 



with the gluon momentum p. This variation can be expressed as [ 17 1 : 



C P (P 2 ) = (-£) 6 I" dtt(l - t) [""J - — , (6) 



37T Jo Jml p — t(l — t)p 



which exhibits the GIM mechanism Wm explicitly. The variation of this 
function with momentum is shown in Fig. 2. Eq. fl5|) is obtained from (|6|) if 
one uses the identification p 2 = — /i 2 , and the approximation m 2 ^> ji 2 3> m 2 . 
For momenta of order Mw, the coefficient Cp is more complicated (T^j, but 



this is irrelevant for us here. Using dimensional regularization Cp(p 2 ) takes 
the following form: 

Cp(p 2 ) = (-^) 6 r(e) r(e) /' dtt(l - t) [(-D u )^ - {-D c )^\ , (7) 

D7T JO 

where T is the gamma-function satisfying T(n) = {n — 1)! for an integer n. 
For space dimension D = 4 — 2e, one has T(e) = 1/e + constant. Moreover, 
r(e) = (— Airfi 2 ) 6 , where /i is a parameter with dimension mass which enters 
the calculation within the dimensional regularization scheme. This parameter 
may be identified with the renormalization scale. The quark masses enters 
the expression (|7]) through the quantity 

D q = m 2 q -t(l-t)p 2 , (8) 

for q = u,c. 

Using the SU(3) matrix relation 



tijtkl = - jj-W , (9) 



splitting the quark current in (^j in a left- and right-handed part and using a 
Fierz transformation, one obtains four different contributions to C(AS = 1). 
The most important penguin four quark operator is Qq, which in the Fierz- 
transformed version is a product of two quark densities, involving both left- 
and right-handed quark fields: 

Q 6 = -8 £K<?i?.) (q R 8 L ) , (10) 

where the sum runs over q = u,d,s. In the naive leading logarithmic approx- 
imation (without RGE), we can do the following identification. 

C 6 = \c P (11) 
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The expression for Cq including RGE can be found in the literature 0, ||, [13 



TJ| . There are also penguin operators which are products of two left-handed 



currents, like 

Qi = 4 ^((I^l) (q L ^s L ) , (12) 
g 

but these are not so important numerically. 

In addition to penguin operators, there is an important pure A J = 1/2 
operator of left-left type which is Q- = Q2 — Qi, where 

Q 1 = 4 (d L ^s L ) (ul1^u l ) and Q 2 = 4 (d L ^u L ) {u L j^s L ) (13) 

are quark operators contained in ([[]). 



3. The Chiral Quark Model 



In the Chiral Quark Model (%QM) JIT], [l2|, [131, chiral-symmetry breaking 
is thought to be taken into account by adding an extra term C x to ordinary 
QCD: 

£qcd — » £qcd x — £qcd + £ x > £qcd = Cf + Cg , 



(14) 



where 



£/ = 9(17 • L> - .M 9 )g ; g 




(15) 



.Mq is the current quark mass-matrix, £ G is the pure gluonic part of QCD, 
and 

C x = -M(qE^qR + m^q L ) ■ (16) 

The constant M in ( ]T6| ) is interpreted as the constituent quark mass, thought 
to be of order 200 - 300 MeV. Note that the constituent and current masses 
are connected to different terms in the Lagrangian, with different transfor- 
mation properties. The quantity X contains the Goldstone-octet fields 7r a : 



Z = exp(iJ2^ a /f) , 



(17) 



where A a are the Gell-Mann matrices, and / = f n — 93 MeV is the pion 
decay constant. The term L x introduces meson-quark couplings. This means 
that the quarks can be integrated out and the coefficients of the various 
terms in the effective meson theory, the chiral Lagrangian, are calculable from 
Cqcdx EH- 1^ has also been found to be suitable for calculating hadronic 



matrix elements of operators obtained from the weak sector [|TJ, |50fl , like in 
eq. ®. 
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There are several versions of the xQM in the literature. Note that in eqs. 
(|14D-(|T7|), there is no kinetic term for the mesons. Thus, the meson fields 
are external. They propagate only after the quarks are integrated out. The 
xQM is thought to apply for momenta of the order and below the scale of 
chiral-symmetry breaking, which we define to be 

A x = 2vr/ 7rV /67iV;, (18) 

where N c = 3 is the number of colours (-numerically A x = 0.83 GeV). There- 
fore a physical ultraviolet cut-off A of the order A x is sometimes used within 
the xQM to parametrize loop integrals which would otherwise be divergent. 
Such divergences will be buried in the physical f v and quark condensate (as 
shown below) , and are not removed by counterterms as in ordinary field theo- 
ries. One should note that A and A x are in principle different quantities. One 
may also use dimensional regularization within the xQM, as shown in [p0| 



Even if dimensional regularization is used, A x is still some kind of effective 
cut-off scale of the xQM, just as for chiral perturbation theory. (Analo- 
gously, the W-mass is some effective cut-off scale of electroweak interactions, 
even if dimensional regularization is used). 

There are a priori several ways to introduce an ultraviolet cut-off A within 
the xQM. The simplest way is to use an explicit (sharp) cut-off in the inte- 
gration over virtual Euclidean momenta. However, this prescription violates 
translation invariance, and will only give a unique result for the leading 



term. One might try a Pauli-Villars type of cut-off [21] , which turns out to 
be technically cumbersome in our case. A good alternative is proper time 
regularization, which has already been used p2[ |15[ in xQM calculations. In 



this case one uses the following replacement for the propagator denominator 
(for Euclidean momenta): 



1 



p 2 + M 2 



/•OO 

/ drexp[-T{p 2 + M 2 )} , (19) 



where £ = 1/A 2 . 

The cut-off A can not be chosen freely. Indeed the x<3^|13|) H3 provides 



a relation between A, M and and eventually gluon condensates. Such 
a relation is obtained because entering the meson- quark coupling ~ 
Mj 5 / f n , is also given by a quark loop diagram for -n — > ^(virtual), and one 
obtains: 

m = ML\f + <^G 2 >+■■■] (20) 

where < ^-G 2 > is the two gluon condensate, and the dots indicate higher 
gluon condensates. The value of A will depend on how many gluon conden- 
sates which are kept in (|20|). In the end both / and will, in the limit 
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iriu,d 0, be identified by f w , but at intermediate stages one might need to 
distinguish them for technical reasons (-we also have fx = f^jp in the limit 
m u ,d, s -> 0). 

The dimensionless quantity f n has the leading behaviour ~ ln(-rW) in a 
cut-off type regularization. Its explicit form will depend on the way the cut- 
off is introduced. Within dimensional regularization, f n = r(e)(47r/i 2 /M 2 ) £ , 
where ft is a parameter with dimension mass, occurring owing to the use of 
dimensional regularization in the \QM. Thus, the would be divergent (for 
A — ► oo) logarithmic term (or 1/e term) is absorbed in the physical pion 
decay constant f w [[ 



The quark condensate is within xQM given by [13, [T5|, |2I], 23 1 



<qq> = -(^/(^^^^ 



N C M 



a 



12M 



< —G > 



TV 



where S(p) is the quark propagator in an external gluon field P3 
tity C q depends on the regularization prescription 



(21) 
The quan- 



-A 2 + M 2 ln AL + 



For a cut-off type reg- 
Within dimensional regularization, 



ularization, C q : 

C q = -M 2 r(-1 + e)(47r/i 2 /M 2 ) e . Equation © puts (in addition to eq.©) 
further restrictions on the parameters of the xQM . 



4. The standard K —> tt amplitude. 



The K —>■ 2n amplitudes can in general be described in terms of chiral 
Lagrangians. Knowing the structure of chiral Lagrangians, the octet part of 
the K — > 27r amplitude can be found from the virtual K —>■ i\ amplitude 
]|, |^, H, [13| . Following |]13[ we define the coupling constant (y^ 1 ^ of the octet 
chiral Lagrangian. In term of this, the virtual K — > n amplitude can be 
written: 

M(K-^n-)s = -V2G F \ u k 2 g ( i /2) f K U, (22) 

where k is the momentum of the virtual meson transition. 

The standard result for the penguin contribution, mainly due to Q 6 , may 



be written as |L3| 



g(Q 6 ) s = g ( i /2) (Q 6 ) s = -16ReC 6 (^) 



(qq) 
f 3 



(23) 



where L 5 is the coupling constant of the relevant term in the strong chiral 
Lagrangian of 0(p 4 ). Within the xQM one obtains WE, EU 



f 3 fi 0) 
iM(qq) 



[1-P] • 



(24) 



7 



where the quantity p is given (up to C(M 4 / A 4 )) by 

a M2 f 1 fort 

P = 6 A^F ; ° r P= f(0^)' (25) 



where the first expression is used within dimensional regularization in [p0|1 , 



and the last one involving the incomplete T-function, T(0,x) with x = 
M 2 /A 2 , is used within proper time regularization in [Tj|. In fact the two ex- 
pressions for p are equivalent when using eq. (|18|), and ( ^0|) with = T(0, x). 
Keeping ffl/f = 1 in the first expression for p in (|25"D, the expression for L 5 
in (|4]) is very sensitive to M. 

Within the xQM, the quantity g£ /2) (Q 6 ) is given by the diagrams in 
Fig. 3. The Ksu and nud vertices are obtained from C x in fll6D-(|T7|). These 
diagrams contain terms independent of the meson momentum. However, 
such terms cancel in accordance with chiral symmetry M . The leading 
physical amplitude corresponds to the sum of the terms of order k 2 . The 
leading term in the expression for L5 in (^) corresponds to diagram 3a, 
while the formally non-leading term ~ p in fl24|) corresponds to diagram 3b. 
Numerically, however, there is a rather strong cancellation between these two 
terms for values of M bigger than 250 MeV. We observe that values of M ~ 
330 MeV or bigger are not acceptable (while keeping fjf'/f = 1) because 
they correspond to zero or even slightly negative values of L 5 in (|24]) . Using 
explicitly the second expression for p in ( p5|) would lead to a reduced value 
of L 5 for M smaller than 250 MeV. One should however take into account 
that r(0, x) = f n is related to the physical f n through (f20|). It is a priori 
an open question if the strong cancellation between the leading diagram 3a 
and the nonleading 3b persists when the momentum variation of the penguin 
coefficient is taken into account (- diagram 3c serves to cancel the momentum 
independent parts of 3a and 3b). 

In |n| the expression ( pi|) was not used. Instead the value L 5 = 1.8 x 10~ 3 
was extracted from other arguments. Moreover, using p ~ M, and the 
scale independent quark condensate < qq >= (— 194 MeU) 3 , it is found that 
g(Qe)s — 0.26. In general, one obtains bigger values for the quark condensate 
if one uses the expression in [|J , 

(mW = — , f ;f* - ( s • (26) 

m u {p) + m d {p) 



Alternatively, using the values|24j] < qq >= (-235MeU) 3 , L 5 = 1.4 x 10 3 



at the scale 1 GeV (- corresponding to M ~ 245 MeV if ([24] ) is used), 



and using C 6 (p = lGeV) = -0.026 for A QC d= 400 MeV, one obtains 



g{Q&)s = 0.15. However, bigger values might be obtained. For p ~ 0.8 GeV 



and Aqcd = 350 MeV, C 6 ~ -0.1, and using in addition p| B§ M ~ 200 
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MeV, one obtains g(Qe)s — 0.8. (Note that, from (|5|) and ( |TT| ) one obtains 
C 6 ~ -0.05 for n ~ A x ) 

The biggest contribution to g^ 2 ^ is coming from the four quark operator 



Q_ of left-left type (-see (|13"D), which gives the amplitude |13| 

\Q-) = lc4i-±r{i-S)} 



(27) 



where C_ is the Wilson coefficient of the operator Q_, which is ~ 2 at \i ~ 
0.8 GeV, and the quantity 5 represents non-factorizable gluon condensate 
corrections: 

N c < ^G 2 > , . 

(28) 
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32tT 2 # 

Writing < — G 2 >= rj 4 , the quantity rj is of the order 400 MeV, which 
gives a value of 5 around 3 (5 ~ 2.6 , 3.0 , and 3.7, for 77 = 376 , 390 , and 
410 MeV, respectively.) Using 77 = 376 MeV and fi ~ 0.8 GeV one obtains 
(Q_) ~ 1.6, while for 77 = 390 MeV and y, ~ 320 MeV one obtains|3| 



9s (Q-) — 2-6. Anyway, the prediction will be below the experimental 
value for the total AI = 1/2 amplitude & IB 



g {1/2) {Tot.) Exp . = 5.1 



(29) 



5. Momentum dependent penguin coefficient 

Now we will present the calculation using a momentum dependent pen- 
guin coefficient. In the chiral limit m Si d — > 0, we obtain the following ampli- 
tude corresponding to the operator Q G from the diagram in Fig. 3a, 

M a (K- - TT; Q 6 ) = v^G^A,, (-8) I x , (30) 

where I x is the two loop integral 

W|f/|f f( ' , ' )[c ^ )1f(,V) - (31) 

Here s and r are loop momenta and fc = r — r' = s — s' is the momentum of 
the virtual K — > 7r transition, and p = s — r is the momentum of the penguin 
gluon. The function F is given by 

r . r > — m 2 

F(r '' r) = p - M2)(r'2 - M2) (32) 
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where the trace in Dirac space has been taken within naive dimensional 
regularization, for simplicity. If Cp is taken as a constant, I x splits in the 
product of two one loop integrals, and we recover the result of other authors[0, 
[|, [L3|, |20fl . Taking Cp as momentum dependent, I x is a two loop integral 
and cut-off procedures are in general cumbersome to implement. We have, 
however, found that if Cp is given as in (^), the integral can be calculated 
analytically within dimensional regularization (in D = 4 — 2e dimensions). 
Then the contribution from the diagram in Fig. 3a given by (|30D , (|3~1~D and (132]) 
turns out to dominate over the contribution from the diagram in Fig. 3b. 

We have calculated the integral using Feynman parametrization in the 
usual way. The r integration is easily performed, while the remaining inte- 
gration over p is the tricky part. Some details will be given in the Appendix. 
Within the approximation m 2 3> M 2 , we find the result: 



c 



Here we have dropped an irrelevant term constant with respect to the virtual 
meson momentum k. This constant cancels with the constant term obtained 
from the diagram in Fig. 3b and 3c. We find that the ~ k 2 contribution from 
diagram 3b is ~ M 2 instead of ~ m 2 . Examining the integral I x to order 
k 2 , one observes (-see Appendix) that one factor 1/e comes from a quadratic 
divergence and will be related to the quark condensate, and the other 1/e 
comes from a logarithmic divergence and will be related to Now we can 
combine the result (|33| ) with the equations (^) and ( ^TJ), and we obtain a 
result for M. a {K~ — > n~;Q G ) in (|30|) . Then we can convert this result into 
the following contribution to g^ 2 ^'- 

^- a fM li+o ^ )+o ^ ] ' (34) 

where the 0(M 2 /A 2 ,) term are coming from 0(e) in (|3^). The result in ([34]) 
has to be compared to the standard amplitude 



3n y A 2 ?>' 2Mf 2 



x J 



(35) 



obtained from combining (p|), (p3|) and (p4|) . It might be surprising that the 
amplitude ([34|) is found to be proportional to m 2 c) which means a huge en- 
hancement compared to (^) from the leading term alone. This behaviour 
reflects that loop momenta of order m c are important in the last loop in- 
tegral over p. The result is very sensitive to the involved parameters. Us- 
ing a s = a s {m 2 ) ~ 0.38 , (qq) = (— 257MeV r ) 3 at /i = m c , we obtain 
98 /2 \Qg) - 10 for M= 0.25 GeV and g(Q 6 ) ~ 19 for M= 0.2 GeV, which 
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is far above the experimental value for g^ 1 ^ 2 ^. However, we observe from our 
calculations that the formally non-leading terms 0(M 2 /A^) in (0), being 
for instance of the form 



M 2 , ml 



. m 2 c 

A 2 <qq> n A 2 



and 



A 2 M 2 

X. 



(36) 



might have sizeable coefficients of order 10. Obviously, big cancellations be- 
tween leading and formally nonleading terms in (£53) has to occur in order to 
get a realistic result. For momentum independent Wilson coefficients, dimen- 
sional regularization works pretty well f20fl . Unfortunately, the coefficients 
(-and their signs!) of the formally nonleading terms in (|36| ) have some ambi- 
guity because dimensional regularization does not distinguish clearly between 
logarithmic divergences (associated with f n ) and quadratic divergences (as- 
sociated with the quark condensate). Because of this, the result (|34D is out of 
control and cannot be used for predicting g{Qo). As a curiosity, we observe 
that while the standard amplitude g(Qe)s has a leading term proportional 
to r(— 1 + e) ■ T(e), the leading term in g(Qe)DR, being the result of a two 
loop integration, is proportional to T(— 1 + 3e) ■ r(— e) ■ F(e), where F(e) is 
finite in the limit e — > 0. 

The quadratic divergence is only appearing in left-right operators like Qq 
in (|H]). From penguin operators containing two left-handed currents, one 
obtains a product of two logarithmic divergences, and there is no ambiguity 
as for quadratic divergences. The standard contribution from Q A is (compare 
with (ff|)) 



A 



(37) 



Using the same method as for the Qq, we obtain for a variable penguin 
coefficient within dimensional regularization (Only Fig. 3a contributes in this 
case, but with axial currents acting at the weak vertices): 



a s / m 2 2 M 2 
^W=— ln^ + --9— 



2 2 

ln^|ln^ + A 

X L X 



(38) 



where A parametrizes some (calculable) non-leading terms, and we have put 



fi = A x in (0). Using the values m c = 1.4 GeV, M ~ 200 MeV and 



A x ~ 830 MeV, we obtain an increase by almost 70% of the leading term with 
respect to fl37j) . However, this is partially compensated by the non-leading 
term ~ M 2 j A 2 X , and the overall result is an increase of order 20 % of g(Qi) dr 
with respect to g(Qi)s- Even if g(Q^) is small (~ 0.05), (j38|) gives an idea 
of how the penguin contributions behave when the variation of the penguin 
coefficients are taken into account. However, g{Q^DR is rather sensitive to 
variations in M. 
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A 


M 


a (A) 


f 


(on) 1/3 


y iV6jp 


y W6jp 


a(Q R ) A r! 
ywej p 


700 


330 


1.00 


101 


— 189 


-0.008 


-0.006 


- 0.01 


830 


330 


0.70 


126 


—222 


0.03 


0.04 


0.03 


1000 


330 


0.53 


155 


-261 


0.10 


0.11 


0.05 


1200 


330 


0.44 


185 


-303 


0.19 


0.22 


0.06 


1400 


330 


0.38 


211 


-342 


0.30 


0.36 


0.07 


700 


250 


1 00 


82 


— 185 


04 


05 


09 


830 


250 


70 


98 


—213 


08 


09 


1 1 


1000 


250 


0.53 


115 


-248 


0.15 


0.16 


0.12 


1200 


250 


0.44 


133 


-285 


0.23 


0.26 


0.13 


1400 


250 


0.38 


148 


-320 


0.30 


0.37 


0.12 


700 


150 


1.00 


47 


-167 


0.09 


0.10 


0.42 


830 


150 


0.70 


53 


-190 


0.11 


0.12 


0.37 


1000 


150 


0.53 


60 


-218 


0.14 


0.16 


0.33 


1200 


150 


0.44 


66 


-248 


0.18 


0.22 


0.31 


1400 


150 


0.38 


72 


-276 


0.21 


0.27 


0.28 



Table 1: Values of y(Q§) P and yi^Q&jp, for constant and momentum depen- 
dent penguin coefficient respectively, calculated for different values of A and 
M. The corresponding modified values g(Q6)pobtained when compensating 
for wrong values obtained for /„- and < qq > are also given. 



Both for Q G and Q 4 , our procedure makes the xQM calculation sensitive 
to m c , which does not sound unreasonable because the scale m c is not too far 
above A x . However, even if matrix elements of the penguin contribution are 
sensitive to the scale m c , we expect that the effect of high momenta ~ m c 
should be damped. Especially, we expect the leading behaviour ~ m 2 c in 
(|34"D to be damped due to some cut-off A. Therefore we have studied the Qq 
amplitude by using proper time regularization. Then the two loop integral 
can not be solved analytically, but we can perform a numerical calculation 
for various values of the cut-off A. This way of regularization will have an 
exponential damping of higher momenta. For this calculation we have used 
© with a s = a, (A). 

The values found for g^ 2 \Q&)p = giQojp for finite cut-offs A calculated 
within proper time regularization are tabulated in table 1. For compari- 
son, the corresponding values g^ 2 \Qe)p = y{Q&)% for a constant penguin 
coefficient (using (||) with fi = A) are also given. As mentioned above, 
the integral for g{Qo) contains the product of a quadratic divergence corre- 
sponding to the quark condensate (see (|2l|)) and a logarithmic divergence 
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for f n (see (pUP). One should note that for various values of A and M, 
the obtained values for f n and < qq > are rather different from the physi- 
cal values (The gluon condensates in ( P0| ) and (^T|) are not included here). 
This leads to wrong values for g(Qo) both for a constant and a momen- 
tum dependent penguin coefficient. To compensate for this, we have calcu- 
lated the modified values g(Qe)p obtained when multiplying with the phys- 
ical /„- and < qq >, and dividing with the respective values f n (A,M) and 
<qq> (A, M) in table 1. This procedure should correspond to the necessity 
of keeping the relation f^> — f as mentioned below (|25|) . For the physical 
quark condensate, we have used < qq > (/i = A) from formula (0), with 
<qq>(fi = lGeV) ~ (— 235MeV A ) 3 , corresponding to (m u + m d ) ~ 12 MeV 
at \i = 1 GeV. 



6. Discussion 

We have investigated the consequences of taking into account the mo- 
mentum dependence of the penguin coefficient for K —>■ 2n decays. The 
calculations are performed within the Chiral Quark Model (xQM). It might 
seem surprising that the penguin amplitude corresponding to Q G gives such a 
big leading contribution within dimensional regularization. But we will point 



out that a similar result was obtained f26[] for the so called Siamese penguin 



diagram, where we found a result ~ m 2 in the limit A 2 — > oo, instead of the 
"expected" ~ A 2 . Unfortunately the dimensional regularization case, which 
is the only one which can be done analytically for Q 6 , is hard to analyze nu- 
merically because dimensional regularization does not distuinguish quadratic 
and logarithmic divergences. One might also worry if dimensional regulariza- 
tion includes too high loop momenta within the xQM. On the other hand, 
the dimensional regularization result for Q 4 seems to be reasonable. 

For finite values of the cut off we have performed a numerical integra- 
tion. Using proper time regularization, we obtain a 10 — 20% increase of the 
penguin contribution compared to the case where the penguin coefficient is 
considered as momentum independent. The reason that we get a slightly 
bigger value for g(Qe)p than g(Qe) P is the following: The relevant diagrams 
for the K — > tt transition to order k 2 are the "eight" diagram (Fig. 3a) and 
the "keyhole" diagram (Fig.3b). For a constant penguin coefficient, it is a 
partial cancellation between these diagrams (in the order k 2 terms) which 
makes the small L 5 in (pi]) . For a momentum dependent penguin coefficient, 
the "eight" and the "keyhole" diagrams are influenced in different ways, and 
g(Qe)p is not anymore directly proportional to L 5 , but to some "smeared 
out" and slightly bigger analogue of L 5 . 

We have not considered the effect of the momentum dependence of the 
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coefficients C± } 2 (or equivalently C±) of the non-penguin operators Qi,2 (Q±) 
in (|13"D . However, these coefficients are not so sensitive to the involved pa- 
rameters, and we expect even smaller effect of momentum dependence than 
for penguin-operators. Furthermore, we have not taken advantage of the 
renormalization group analysis involving the penguin coefficient Cq. This is 
because the individual quark momenta can hardly be distinguished, owing 
to the mixing of operators. 

Our analysis offers no solution to the AI = 1/2 puzzle for K —* 2tt 
decays. Moreover, the 10-20% increase of the penguin contribution seen in 
table 1 is moderate compared to the increase in the penguin coefficient C§ 



owing to two loop effects |19| . The penguin contribution to K — > 2tt is very 
sensitive to the involved parameters as the quark condensate, the constituent 
quark mass M, and the coefficient Cq through Aqcd and the matching scale 
jj,. In conclusion, other contributions than (Qe) must be important for 



the AI = 1/2 enhancement |25| 



7. Appendix 

We will here give some details from the calculation of I x in (|31~|) . We 
expand the propagators up to second order in the external virtual meson 
momentum k. Then the propagators [r 2 — M 2 ] -1 and [(r — p) 2 — M 2 ]^ 1 
occur in some powers. We need one Feynman parameter (x) to perform the 
integral. We use the general formula 



(-tt) d / 2 T(n-D/2) 



where 



Jn(D) 



(39) 



(40) 



{2ir) D T{n) 

Note that J n+ i(D) = J n (D)(n — D/2)/n. Dropping the constant term ~ k° 
we obtain 



iJ 3 (D)k 2 J o dx (s c (3) 



6 
D 

(£^£)(l-2x(l-x)) 

(c -»• u) , 



l-x(l-x)) -2 



,Q-D 
D 



)x{l-x) M 2 S C (A) 



(41) 



where the quantity S q (n) depends on the current quark masses for q = c,u: 



S q (n) 



d D p 
(2^p 



[Cp(p 2 )] q [-A] 



(D/2-n) 



(42) 
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where [C P (p 2 )] = [C P (p 2 )} u - [C P (p 2 )} c , and A = M 2 - x (1 - x) p 2 . To 
find S q (n), we have to use the dimensional regularization version of [Cp{p 2 )} 
in (0), and we need the generalized Feynman parametrization 

jtttt; (43) 



(-A)'(-£>«)« B(i,c) Jo Nj+ 

where I = n - D/2 and iV g = (1 - y)(-A) + y(-D q ) = h(p 2 - WT 2 ). Here, 
/i = yt(l —t) + (l —y)x(l —x) and m^ 2 = (ym 2 + (1 —y)M 2 )/h. The quantity 
-B is given by 



S(a,6) = / das " 1 (1- re) 6 " 1 = 7 ^ (44) 



r(q) ■ m 

o "~~ ' ' ' Via ■ b) 

Using again d39|) , we obtain 

r ( e ) /-S^/i ,n Z 1 j (i-y)'-y- 1 



S q (n) = iJ^D)^Lf o dtQt(l-t)J o dy [ - {-m?Y , 

(45) 

where u = D/2 — (/ + e). Now, for the charm quark case g = c, we use the 
approximation m^T 2 = ym 2 /h. Then we obtain 



5 c (n) = iJ l+e (D)^P^{-m 2 c ) u f dt6t(l - t) I{x,t) (46) 
£>{l,e) Jo 



Then the clue is to recognize the integral 

y) l - l y e+u - 1 

h u + l+ * " [x(l - x)] l [t(l - t)} 



which gives 



S c (n) = 3 c (n)[x(l-x)]- 1 (48) 

with 

r(c) 



5 c (7i) = iJ l+e (D)-^-(-m 2 c ) u B(l,u + e)6B(2-u-e,2-u-e) (49) 
Then, to leading order in m 2 , we find 



I x = -iJz(D)k 2 S c (?>)B(D/2-2,D/2-2) (50) 

Here the factor B(D/2—2, D/2— 2) corresponds to the logarithmic divergence 
connected to /„-, and Ji +e (D) in S c (3) to the quadratic divergence connected 
to the quark condensate. 
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Figure 1: The penguin loop diagram. At the lower line q — u,d, s. 



19 



P [p*2] 
10 r 



8 



6 



2 ■ 




0.5 ' ' ' ' ' 1.5 ' ' 2 ' ' ' 2 .' 5 ' ' ' 3 P 



Figure 2: Variation of the penguin coefficient Cp = —^P(p 2 ) in (6) with 
squared Euclidean momentum p 2 (-in units GeV 2 ). The horisontal line rep- 
resents the value in eq. (|J) for /i = 0.83 GeV. 
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Figure 3: Diagrams for the penguin operator Q 6 contribution to K — > ir. 
The solid lines represent quarks. Diagram (a), upper left, is the "eight" dia- 
gram, while (b) is the "keyhole" diagram. Diagram (c) serves to cancel the 
momentum independent part of (a) and (b). 
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